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A spin- and time-dependent electron transport has been studied in a paramagnetic resonant 
tunneling diode using the self-consistent Wigner-Poisson method. Based on the calculated current- 
voltage characteristics in an external magnetic field we have demonstrated that under a constant 
bias both the spin-up and spin-down current components exhibit the THz oscillations in two different 
bias voltage regimes. We have shown that the oscillations of the spin-up (down) polarized current 
result from the coupling between the two resonance states: one localized in the triangular quantum 
well created in the emitter region and the second localized in the main quantum well. We have also 
elaborated the one-electron model of the current oscillations, which confirms the results obtained 
with the Wigner-Poisson method. The spin current oscillations can lower the effectiveness of spin 
filters based on the paramagnetic resonant tunneling structures and can be used to design the 
generators of the spin polarized current THz oscillations that can operate under the steady bias and 
constant mae: netic fieldQ 
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I. INTRODUCTION 

Dilute magnetic semiconductors (DMSs) such as Zn- 
MnSe and GaMnAs are promising materials due to 
their potential applications in spintronics.'^l'^ The re- 
cent homo- and heteroepitaxy method^^HI^ allow to de- 
posit the DMS layers with thicknesses of few nanometers, 
which enables to fabricate spintronic nanodevices with 
the spin polarization of the current being controlled by 
the external magnetic and electric fields. 

The spin filter based on the DMS was proposed by 
Euges.!^ The nanostructure studied in Ref. il6j con- 
sisted of a paramagnetic semiconductor layer made from 
Zni-a^Mn^^Se sandwiched between two non-magnetic 
ZnSe layers. The external magnetic field leads to the 
giant Zeeman splitting of the conduction band minima 
in the paramagnetic layer, which causes that this layer 
acts as a potential well for spin-down electrons and a 
potential barrier for spin- up electrons. As a result, the 
total current flowing through the nanostructure is domi- 
nated by the spin-down electrons. In such a spin filter, 
the change of the spin polarization of the current re- 
quires the change of the external magnetic field. In 
the recently fabricated paramagnetic resonant tunneling 
diodes (RTDs]^^, the electrical control of the spin po- 
larization of the current has been achieved. The spin- 
polarized current is controlled by the bias voltage in the 
prese nce of the external magnetic field in paramagnetic 
j^rpj-^ji7|i8|20|2T| even without the external magnetic 
field in ferromagnetic RTDs.^SlHUl j£ ^j^^ quantum well 
in the RTD is made from the DMS, the spin splitting 
of the quasi-bound state energy level gives rise to the 
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resonant tunneling conditions for the spin-up and spin- 
down electrons satisfied for different bias voltages. This 
leads to the separation of both the spin current compo- 
nents and consequently to the spin polarization of the 
net current. The operation of the paramagnetic RTD 
based on ZnSe/ZnBeSe/ZnMnSe heterostructure as a 
spin filter has been experimentally demonstrated by Slo- 
bodskyy et al.^^ and theoretically described by Havu et 
al.^^ Recently, the spin polarization of the current in the 
paramagnetic RTD has been reported at zero magnetic 
field. All these studies of the spin-polarized currents 
in the magnetic RTDs were devoted to the stationary 
(steady) currents. The oscillations of the spin-polarized 
currents in magnetic RTDs have not been studied un- 
til now, although the intrinsic oscillations of the current 
have been detected in the non-magnetic resonant tunnel- 
ing structures .I^^M^ 

The intrinsic current oscillations occurring in the non- 
magnetic RTDs are intensively studied due to their po- 
tential application as THz generators. Recently, the THz 
oscillations have been experimentally observed in the 
Gain As/ Al As RTD integrated with a slot antenna.'^^M^ 
The attempts undertaken in order to explain the origin of 
these oscillations led to contradictory conclusions. Ricco 
and Azbel^^ argued that the oscillations of the current in 
the non-magnetic RTD result from the fact that the sys- 
tem enters and leaves the resonant current conditions. 
According to this model, the intrinsic oscillations of 
the current should occur at the resonance bias that is in 
contradiction with the experimentaP^M^and numericaP^ 
studies, which show that the current oscillations occur 
only in the negative differential resistance (NDR) regime 
of the cur rent- volt age characteristics, i.e., above the res- 
onance bias. Woolard et al.^^ suggested that the cur- 
rent oscillations in the non-magnetic RTD result from the 
charge fluctuations in the potential well created between 
the emitter and the nearest barrier. This proposition 
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was extended by Zhao et al.p^M^who showed that the 
intrinsic oscihations are due to the couphng between the 
quasi-bound state locaUzed in the emitter-related quan- 
tum weh and the quasi-bound state in the main quantum 
weh. However, this model^^ is based on the adiabatic 
approximation, according to which the electron states are 
slowly varying in time. This assumption is not valid in 
the THz oscillation regime. Moreover, the explanation 
given by Zhao et al. l^^^ does not answer the question, 
why the current oscillations do not decay in time as a re- 
sult of the dissipative factor corresponding to the imag- 
inary part of the resonance state energy. In our recent 
paper, we have shown that the intrinsic oscillations of 
the spin-unpolarized current result from the coupling be- 
tween the two arbitrary quasi-bound states in the non- 
magnetic triple-barrier resonant tunneling structure. We 
have found the intrinsic oscillations in two bias ranges. 
We have interpreted the oscillations that occur at the 
bias values below the resonance bias as resulting from 
the coupling between the quasi-bound states in both the 
quantum wells, while the oscillations that occur in the 
NDR regime as resulting from the coupling between the 
emitter-related quasi-bound state and the quasi-bound 
state localized in the nearby quantum well. 

In the present paper, we study the spin- and time- 
dependent electronic transport in the paramagnetic 
double-barrier RTD based on ZnSe/ZnBeSe/ZnMnSe. 
We show that - under certain well-defined constant bias 
and external magnetic field - the intrinsic oscillations of 
both the spin current components appear, which leads 
to the oscillations of the spin polarization of the total 
current. Based on the analysis of the spin- and time- 
dependent potential energy profiles and electron density 
distributions, we demonstrate that the oscillations of the 
spin polarized current result from the coupling between 
the quasi-bound state localized in the triangular quan- 
tum well created in the emitter (E) region with the spin- 
dependent quasi-bound state localized in the main quan- 
tum well (QW). The present results show that at certain 
bias voltages the spin polarization of the current fiow- 
ing through the paramagnetic RTD is not constant but 
oscillates with the THz frequency. This effect can be of 
crucial importance for the possible applications of this 
nanodevice as a spin filter and THz generator. 

The paper is organized as follows: in Sec. II, we de- 
scribe the theoretical model of the paramagnetic RTD 
and the time-dependent Wigner-Poisson method. Sec- 
tion HI contains the results. Section IV - the discussion, 
and Section V - conclusions and summary. 



II. THEORETICAL MODEL 

We consider the spin-dependent electronic transport 
through the paramagnetic RTD that consists of the para- 
magnetic quantum- well layer made from Zni-a^Mn^^Se 
embedded between the two barrier layers made from 
Zno.95Beo.05Se (Fig. [1]). The active (undoped) region 



of the nanodevice is separated from the n-doped ZnSe 
ohmic contacts by the two spacer layers made from 
ZnSe. In the presence of the external magnetic field 
B = (0, 0, B) applied in the growth (z) direction, the 
exchange interaction between the conduction band elec- 
trons and the Mn^+ ions leads to the giant Zeeman split- 
ting of the conduction band minimum in the paramag- 
netic quantum- well layer. ^ In the external magnetic field, 
the conduction-band electrons form the Landau states 
with the wave functions spread over the x — y plane. 
Therefore, the one-electron problem can be separated 
into (x, y) and z coordinates and the electronic trans- 
port between the emitter and collector can be described 
as the one-dimensional motion of the electron in the z 
direction. 




FIG. 1. (Color online) Potential energy in the paramagnetic 
RTD for the electrons with spin- up (solid, red) and spin-down 
(dotted, blue). Coordinate z is measured along the growth 
direction, iie{C) is the electrochemical potential of the emit- 
ter (collector). The active (undoped) region consists of the 
paramagnetic quantum well made from Zno.92Mno.08Se sand- 
wiched between the two Zno.95Beo.05Se potential barriers and 
is separated from the n-doped ZnSe ohmic contacts by the 
two ZnSe spacer layers. 

In order to simulate the electronic transport through 
the nanostructure, we apply the time-dependent Wigner- 
Poisson approach, according to which the conduction 
band electrons are described by the spin-dependent 
Wigner distribution function.'^^^^^ The quantum kinetic 
equation takes on the following fornP^ 



dp^{z,k,t) , nkdp^{z,k,t) 



dt 



m 



dz 



(1) 



2^h 



+00 



.k-k',t)p^{z,k',t) , 



where p^{z^k,t) is the spin-dependent Wigner distribu- 
tion function, k is the z-component of the wave vector, 
m is the electron conduction-band effective mass, and 
(J = (t, i) is the spin index. 
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The non-local potential hicj{z^ k — k']t) for spin channel 
<j is given by the formula 



+ CX) 



Ua{z,k-k',t)= dz' [Uaiz + z' /2,t) 



(2) 



-Ua{z- z'/2, t)] exp [-i{k- k')z'] , 

where Ua{z,t) is the spin-dependent potential energy, 
which can be expressed as the sum of the two terms 



U,{z,t) = U^\z,B) + U,i{z,t). 



(3) 



In Eq. (|3|, the first term denotes the spin-dependent po- 
tential energy of the conduction-band bottom, while the 
second term is the potential energy of an electron in the 
electric field acting in the nanostructure and has the form 



Uel{z,t) = Uv,{z)^UH{z,t) 



(4) 



where Uvi,{z) is the potential energy of an electron in the 
electric field generated by voltage applied between the 
emitter and collector and UH{z;t) is the Hartree energy 
that takes into account the electron-electron interactions. 
In Eq. (|3|, we have neglected the exchange energy (see 
Subsection IV.C). 

The conduction-band potential energy profile has the 
form 

/7o , zi < z < Z2 and zs < z < Z4, 

U^\z, B) = { ±Ez{B) , ifz2<z< zs, 







otherwise, 



(5) 

where Zi and Z2 {zs and Z4) are the positions of the left 
(right) potential barrier interfaces, Uq is the height of 
the potential barrier, and Ez{B) is the Zeeman energy 
of an electron in the paramagnetic quantum- well layer, 
which is positive (negative) for the spin-up (spin-down) 
electrons. For a small concentration of Mn^+ ions the 
Zeeman energy can be expressed by the formula^ 



Ez{B) 



2 V ksTeff J 



(6) 



where Noa = 0.26 eV is the sp-d exchange constant, x 
is the concentration of Mn^+ ions, Bs is the Brillouin 
function for spin 5 = 5/2 that corresponds to the spin 
of Mn^+ ion, g is the effective Lande factor, /i^ is the 
Bohr magneton, 5*0 and Tg// are the phenomenological 
parameters corresponding to the antiferromagnetic inter- 
action between the Mn^+ ions and have been taken on 
as S'o = 1.18 and Teff = 2.55 K for x = 0.083."^^ 

Potential energy Uei{z, t) satisfies the Poisson equation 



dz2 



-[ND{z)-n{z,t)] , 



(7) 



where e is the elementary charge, Eq is the vacuum elec- 
tric permittivity, e is the relative static electric permit- 
tivity, Nd{z) is the concentration of the ionized donors. 



and n{z) = J2a'^(^(^) total density of electrons 

with ncj{z) being the density of the electrons with spin 
(J. The energy of the emitter conduction-band bottom is 
taken as the reference energy and set equal to 0. 

In order to solve the system of non-linear equations 
([1]) and ([jl) we apply the self-consistent time-dependent 
procedure^^ with the boundary conditions Uei{0,t) = 
and Uei{L,t) = —eVb for the Poisson equation, where L 
is the length of the nanodevice. For the quantum ki- 
netic equation ([T]) we apply the generalized form of the 
boundary conditions proposed by Frensle^BS 



p^{L,k,t) 



(8) 



k>0 



k<0 



Distribution function f^{k) {v = E,C) is derived from 
the Fermi-Dirac distribution function by summing over 
the Landau energy levels, which leads to 



eB v-^ 



n=o exp 



1 



\ 2m ^na J 



(9) 



where T is the temperature, E^^ = jUjy — hujc (n + 1/2) — 
(jIIbB, fijj = Fj, — eVu is the electrochemical potential 
of reservoir Fy is the corresponding Fermi energy, Vy 
is the voltage applied to contact and ujc = eB/m is 
the cyclotron frequency. In Eq. ([9|, N^y^ax determines the 
highest occupied Landau state and is determined by the 
density of electrons in the reservoirs. 

The spin-dependent Wigner distribution function al- 
lows us to calculate the spin-dependent electron density 



+00 



{z,k,t) 



(10) 



and the spin-dependent current density 



2ttL 



L +00 



m 



(11) 



The spin polarization of the current is defined as follows: 



Pit) 



Mt)-jdt) 



(12) 



where and are the spin-up and spin-down current 
densities, respectively. 

In the present work, we focus on the oscillations of the 
spin-polarized currents. In order to choose the nanos- 
tructure parameters, for which the current oscillations 
are mostly pronounced, we have performed a number 
of computational runs that allowed us to establish the 
optimal values of the parameters. The present simula- 
tions have been performed for the following values of the 
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FIG. 2. (Color online) Current- voltage characteristics for spin-up (curves with circles, red) and spin-down (curves with triangles, 
blue) electrons and for the magnetic field (a) B = 2 T, (b) B = 4 T, and (c) B = 6 T. The hatched areas correspond to the bias 
regimes, in which the oscillations of the spin-polarized current occur. The values of the current density in these regimes are 
calculated by averaging over one period of oscillations. Inset (A) shows how the steady current is reached after several initial 
fluctuations for the bias denoted by arrow (A) on main panel (a), inset (B) shows the formation of the stable oscillations of 
the spin currents that occurs in the bias regimes denoted by arrows (B) on main panel (a). 



nanostructure parameters: the thickness of each contact 
is equal to 17 nm, the thickness of each spacer layer is 
3 nm, the thickness of the potential well layer is 5 nm, the 
total length of the nanodevice L = 54 nm, and the Mn 
concentration x = 8.3 %. The barriers are assumed to be 
symmetric with thickness 3 nm each. We have taken on 
the height of both the potential barriers Uq = 0.115 eV.^^ 
The contacts are made from the n-type ZnSe with the 
homogeneous concentration of donors Njj = 10^^ cm~^. 
Due to the small thickness of the double-barrier region 
we assume the effective electron mass of the ZnSe con- 
duction band, i.e., m = 0.16 mo, where mo is the free 
electron mass. We take on the relative electric permit- 
tivity of ZnSe e = 8.6 for the entire nanostructure. The 
present simulations have been carried out for tempera- 
ture T = 1.2 K on the computational grid with = 95 
mesh points for coordinate z and = 72 mesh points 
for z-component of the wave vector. The z coordinate in- 
crement was Az = a, where a = 0.5667 nm is the lattice 
constant of ZnSe. 



III. RESULTS 

The calculated current- voltage characteristics for spin- 
up and spin-down components of the current in the pres- 
ence of the external magnetic field are displayed on Fig.[2j 
The resonant current peaks for both the spin current 
components are separated, which results from the fact 
that the resonance conditions for the spin-up and spin- 
down electrons are satisfied at different bias voltages due 
to the giant Zeeman splitting of the quasi-bound state 
energy levels in the paramagnetic quantum well. If the 
magnetic field increases the resonant current peak for the 
spin-down electrons shifts towards the lower bias, while 
the resonant peak for the spin-up electrons shifts towards 
the higher bias [Fig. [2|a,b,c)]. The separation of these 
peaks is the basic property that can be exploited to fab- 
ricate an effective spin filter, in which the spin polar- 
ization of the current is controlled by the bias voltage. 
Figure [3] shows the spin polarization of the current as a 
function of the bias voltage for different magnetic fields. 
For all values of the magnetic field the spin polarization 
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of the current varies from P = — 1 for the low bias to 
P :^ +1 for the sufficiently high bias. Moreover, for the 
higher magnetic field the transition between the almost 
fully spin polarized currents occurs in a narrower bias 
interval. 



Q- 0.0 




0.12 



FIG. 3. (Color online) Spin polarization P of the cur- 
rent as a function of bias voltage Vh for the magnetic fields 
5 = 2, 4, 6 T. For the oscillating current the spin polariza- 
tion has been calculated by averaging over one period of the 
oscillations. 

The current- volt age characteristics presented in Fig. |2] 
have been calculated in the bias range from V5 = to 
Vb = 0.2 V with step dVt = 0.005 V. For each step, 
the bias voltage has been changed to the next value by 
SVb if both the spin current components reach the steady 
state as shown in inset (A) of Fig. [2|a). We have found 
that at certain bias voltages that belong to the NDR 
regimes for the spin-up and spin-down electron currents 
both the spin currents do not stabilize as the steady cur- 
rents but oscillate with the THz frequency [cf. inset (B) 
on Fig. Ha)]. 

In Fig. [4j we present the formation of the current os- 
cillations for the bias Vb = 0.092 V that corresponds to 
the NDR regime for the spin-down current component 
[Fig.[4];a)] and Vb = 0.112 V that corresponds to the NDR 
regime for the spin-up current component [Fig. |4]^b)]. 
In both the cases, the spin-polarized current oscillates 
with constant frequencies equal to fosc = 3.7 THz for 
Vb = 0.092 V and fosc = 3.2 THz for V5 = 0.112 V. 
For the bias voltages, for which the current oscillations 
have been found, the spin current polarization is not well 
defined but also oscillates with the THz frequency. Fig- 
ure [5] displays the oscillations of the spin current polar- 
ization that occur in the NDR regime for the spin-down 
[Fig.jsja)] and spin-up [Fig.jsjb)] current component. We 
note that the spin current polarization exhibits the oscil- 
latory behavior in the bias regimes, in which it reaches 
the maximal values. 

In order to get a physical insight into the mechanism of 
oscillations of the spin-polarized currents, we have calcu- 
lated the time- and space-dependence of the electron den- 
sity in the different regions of the nanostructure (Figs. |6] 



■ (a) V = 0.092 V, B = 2 T, f =3.7 THz 

\ / b ' 'osc 

c 0.10 - Ji'.jiJuV'* 



-I — ' — I — ' — I — y/^-r- 




10 20 30 58.2 58.8 59.4 60.0 
time (ps) 



0.16 
0.14 
0.12 

0.10; 

0.03 
0.02 
0.01 
0.00 




(b) V^= 0.1 12 V, B = 2 T, f^= 3.2 THz 



- spin-up 
• spin-down 




4 6 8 10 
time (ps) 



12 14 



FIG. 4. (Color online) Spin-up (solid, red) and spin-down 
(dashed, blue) current oscillations in the NDR regimes for 

(a) the spin-down current component (Vb — 0.092 V) and 

(b) for the spin-up current component (Vb = 0.112 V). The 
calculations have been performed for B — 2T. 



and[7|). Figure [g] shows the time dependence of spin-up 
and spin-down electron densities in the emitter region 
and in the main quantum well (QW) calculated for the 
bias corresponding to the NDR regimes for spin-down 
and spin-up current components. We see that both the 
spin-up and spin-down electron densities oscillate with 
the constant frequencies that are equal to the frequen- 
cies of the oscillations of the corresponding spin-polarized 
current. Figure |6|a) shows that in the NDR regime for 
the spin-down current component, the amplitude of the 
spin-down electron density oscillations in the QW region 
is much larger than the corresponding amplitude of the 
spin- up electron density oscillations. At time instant ti, 
the spin-down electron density reaches the maximal value 
in the QW region and the minimal value in the emitter 
region. On the other hand, at time instant t2, the spin- 
down electron density is minimal in the QW and maximal 
in the emitter region [cf. Fig. |6ja)]. The similar behav- 
ior has been obtained for the spin-up electron density in 
the corresponding NDR regimes for the spin-up current 
component [Fig. |6jb)]. 

Figure [T] depicts the spatial distribution of the spin- up 
and spin-down electron densities together with the cor- 
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(a) V^=0.092 V, B = 2 T, f„^=3.7 THz 
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FIG. 5. (Color online) Oscillations of the spin current po- 
larization P in the NDR regime for (a) the spin-down current 
component {Vh — 0.092 V) and (b) spin-up current compo- 
nent (Vb = 0.112 V). Insets display the regions marked by the 
rectangles on the main panels. 



responding self-consistent potential energy profiles deter- 
mined at time instants ti and t2 marked on Fig.[6j Figure 
[7|^a) shows that the density of the spin-down electrons 
in the QW is fairly large at time instant ti and much 
smaller at ^2- These changes of the spin-down electron 
density in the QW mean that the tunneling conditions at 
time instants ti and t2 are different. Therefore, the spin- 
down current component approaches (at ti) and partially 
leaves (at ^2) the resonant tunneling conditions. At these 
time instants, the density of the spin-up electrons in the 
QW is considerably smaller and oscillates with smaller 
amplitude [Fig.[7|^a)]. In this case, the low-amplitude os- 
cillations of the spin-up electron density are induced by 
the high-amplitude oscillations of the spin-down electron 
density [cf. Fig. (joja)]. Similar results have been ob- 
tained for the spin-up electron density in the QW for the 
bias corresponding to NDR regime of the spin-up current 
component [Fig. [7jb)]. At time instant ti the spin-up 
electron density in the QW is larger than at time instant 
t2. At these time instants, the density of the spin-down 
electrons in the QW is considerably smaller and changes 
only slightly during the oscillations [Fig. [7|^b)]. This in- 
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FIG. 6. (Color online) Oscillations of the spin-up (red) 
and spin-down (blue) density of the electrons localized in the 
emitter region (dashed line, right scale) and the QW (solid 
line, left scale) calculated for magnetic field 5 = 2 T and the 
bias voltage from the NDR regimes for (a) spin-down current 
component (Vb = 0.092 V) and (b) spin-up current component 
(Vb = 0.112 V). The density of electrons localized in the QW 
reaches the maximum (minimum) at time instant ti (t2). 



dicates that the spin-down polarized current is out of the 
resonance, which explains the small values of the spin- 
down current in the considered bias regime. 



IV. DISCUSSION 



The results of Sec. Ill show that - for the paramag- 
netic RTD in the external magnetic field - we can ex- 
pect two bias voltage regimes, in which the spin polar- 
ization of the current is not well defined but oscillates 
with the THz frequency. Figure [2] shows that the os- 
cillations of the spin-polarized current are generated in 
the two following bias ranges: the first corresponding to 
the NDR regime for the spin-down current component 
and the second corresponding to the NDR regime for 
the spin- up current component. The present results ob- 
tained for the paramagnetic RTD resemble those for the 
non-magnetic RTD,^^ according to which the current os- 
cillations occur only in the NDR regime of th e cur rent- 
voltage characteristics. The previous 
the current oscillations in the non-magnetic RTD allow 
us to conclude that the occurrence of this phenomenon 
does not depend on the spins of electrons flowing through 
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FIG. 7. (Color online) Density of the spin-up (red) and 
spin-down (blue) electrons and the corresponding profiles of 
the self-consistent potential energy as a function of coordinate 
z at time instants t\ and marked on Fig. |6] for bias (a) 
H = 0.092 V and (b) Vb = 0.112 V. 



the nanostructure. In the next subsections, for simplic- 
ity and without loss of generality, first we propose the 
model of the spin-independent current oscillations in the 
non-magnetic RTD and next extend this model to the 
description of the spin-polarized current oscillations in 
the paramagnetic RTD. 



A. Model of current oscillations 

In this subsection, we will analyze the current oscilla- 
tions using the one-electron approximation. In order to 
explain the origin of the oscillations, let us consider the 
physical processes occurring in the NDR regime. In this 
regime, the system leaves the resonant tunneling condi- 
tions, which means that the probability of the reflection 
of the electron from the emitter barrier becomes fairly 
large. The interference between incident and reflected 
electron wave functions together with the bias-induced 
lowering of the electron potential energy leads to the 
formation of a shallow triangular quantum well in the 
emitter region (cf. Fig. [t]), in which the quasi-bound 
(resonance) state is formed. Let us consider the coupling 
between the resonance state = i\)e mainly localized 
in the triangular quantum well in the emitter region and 
the resonance state ?/^2 = i^QW localized in the QW. In 



the framework of the proposed model, we will calculate 
the probability of transition from state to state i\)2 
under influence of the weak perturbation corresponding 
to the change of the bias and show that this probability 
oscillates in time. 

The one-electron wave function ^(z^t) evolves in time 
according to the time-dependent Schrodinger equation 



2m dz^ 



^(z,t) + /7[^(^,t)]^(z,t), (13) 



where potential energy [/[^(z,t)] is a functional of the 
wave function, which results from the wave-function rep- 
resentation of the electron density n = n[?/^(z,t)] that in 
turn determines the potential energy via Poisson equa- 
tion ([7|. Potential energy U[^{z^t)] can be expressed as 
the sum of three terms 

U[^{z,t)] = U\z)^SUv,{z)^SUh[^{zM (14) 



where U^{z) is the potential energy of the electron in 
the steady state for the bias = V^^^ that corresponds 
to the maximum of the resonant tunneling current. The 
Hamiltonian with potential energy U^{z) describes the 
unperturbed system. The next two terms in Eq. (14) cor- 
respond to the perturbation, where SUvh{z) is the change 
of the electron potential energy caused by varying the 
bias from V5 = V;,^"" to V5 = V^"-^' + SV^ {6Vb > 0) and 
dUH[^{z^t)] is the corresponding change of the Hartree 
potential energy. 

The change of Hartree potential energy SUH[^{z^t)] 
results from changing the potential energy of the elec- 
tron in the external electric field by amount SUvk{z) and 
occurs with some time delay with respect to dUvi,{z). 
Moreover, SUH[^{z,t)] is much smaller than 6Uvi,{z), 
and therefore, in the first approximation, it can be ne- 
glected when considering the electron states at the bias 
Vb = V;-'+6Vb. 

In order to avoid the exponential divergence of the res- 
onance wave function in Eq. (13) all the calculations will 



be performed using the complex-scaling theory of reso- 
nance states. The complex-scaled Hamiltonian of the 
electron has the form 



H{ze'^) = 



2m 



-^U''{ze'')^SUv,{ze'') = 

(15) 



dz^ 



where 6 is the scaling parametei^^ and i^o(^e*^) is the 
Hamiltonian of the unperturbed system for the bias V5 = 



T/res 



In the NDR regime, i.e., for = V^^^ + (514, the reso- 
nance states '0i(^) and V^fl^) almost degenerate, i.e.. 
El = E2 -\- SE. Therefore, the one-electron wave func- 
tion can be expressed as a linear combination of the wave 
function '^1(2;) being localized in the emitter region and 
V^2(^) localized in the QW 



^^(z,t) =ai(t)V^f(z)+a2(t)V^2^(z). 



(16) 
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Wave functions '^^(2;) (n = 1,2) are calculated using 
the time-independent Schrodinger equation for the un- 
perturbed Hamiltonian 



(17) 



where the energy of the resonance state En = — iTnf^ 
consists of the real part that determines the position 
of the resonance on the energy scale and the imaginary 
part that defines the width of the resonance and de- 
termines its lifetime. 

After inserting (|T6| into the time-dependent 
Schrodinger equation (|13|) we obtain 



/ daijt) 
\ dt 



where 



El + Wii W12 
W21 E2 + W22 



ai{t) 
a2{t) 



(18) 
(19) 



Finally, the wave function takes on the form 

ii,0{z,t) = (a;e-t«+* + aOe-*«-*) <(^) 

+ (a?a+e-^«+* + a^a-e-*«-*j ^/;^(z), (20) 
where 2 are the time-independent coefficients, 

W21 



a 



u - W22 



■ E2 + Wii + W22 T V A 



and 



A=iEi+ Wu -E2- W22f 



(21) 



(22) 



(23) 



If we take the resonance state i^i{z) as the initial state 
of the system, i.e., the state for t = to and V5 = V^^^ ^ 
the probability of transition from state il^i{z) localized in 
the emitter quantum well to state V^fC^) localized in the 
QW is expressed by the formula 



^ sin^ u)t 



where 



w = 



2n 



(24) 



(25) 



The probability of the electron transition from the QW 
to emitter is given by 



2^1 



it) = 1 



1^2 



{t). 



(26) 



Based on Eqs. (24) and (26) we state that the coupling 



in the QW (measured by W12) causes that the probabili- 
ties of the electron transition between the both resonant 
states oscillate with frequency uj. As a consequence of 
this coupling, both the electron density and the current 
density oscillate in certain well-defined bias intervals in 
the NDR regime. Eqs. ([24]) and ([26]) show that the den- 



between the resonance states localized in the emitter and 



sities of electrons localized in the emitter and in the QW 
regions oscillate in anti-phase, which supports the results 
presented in Fig. [6] 



Spin current oscillations in the paramagnetic 
RTD 



Based on the results of Sec. Ill and Subsection 
IV. A we argue that the oscillations of the spin-polarized 
currents in the NDR regimes result from the coupling 
between the spin-dependent quasi-bound state local- 
ized in the emitter region with the corresponding spin- 
dependent quasi-bound state localized in the main QW. 

Now we will analyze the formation of the oscillations of 
the spin-down polarized current in the NRD regime of the 
current-voltage characteristics for this current compo- 
nent [cf. Fig. [2]^a)]. In this regime, the resonant tunnel- 
ing conditions for the spin-down electrons are no longer 
satisfied and the refiection probability of these electrons 
from the emitter barrier becomes fairly large. Due to 
the interference and the bias-induced lowering of the po- 
tential energy, the spin-down quasi-bound state of the 
electron is formed in the triangular quantum well in the 
emitter region close to the left barrier [cf. Fig. [7]^a)]. 
We have estimated the energy of this resonance state to 
be E'^ —0.1 meV. This energy differs by few meV 
from energy Eq^^ of the resonance state localized in 
the QW. Therefore, these states can be treated as al- 
most degenerate and we can apply the results of Sub- 
section IV. A. The coupling between these two quasi- 
bound states leads to the oscillatory transitions of the 
spin-down electrons between the emitter and QW regions 
[cf. Fig.[6];a)] described by the Eqs. ([24| and (|26|. The 
oscillatory changes of the spin-down electron density in 
the QW causes that the bottom of the QW oscillates 
in time, which gives rise to the oscillatory changes of the 
tunneling conditions. Figure [sj^a) shows the transmission 
coefficient T as a function of incident electron energy E 
at time instants ti and ^2 marked on Fig. [6]^a). The po- 
sitions of the peaks of T{E) correspond to the energies 
of the spin-up and spin-down resonance states localized 
in the main QW. Figure [8]^a) shows that at time instant 
ti the resonance tunneling condition is satisfied for the 
spin-down electrons, which leads to the accumulation of 
these electrons in the QW [Fig. [7]^a)]. After ~1 ps the 
system gets out of the resonance and the density of the 
spin-down electrons in the QW decreases reaching the 
minimal value at time instant ^2 [cf. Fig. [7]^b)]. As a 
consequence, the spin-down polarized current oscillates 
[cf. Fig. ([6](a)]. For this bias voltage, the oscillations of 
the spin-up electron density are induced by the electron- 
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FIG. 8. (Color online) Transmission coefficient T as a func- 
tion of incident electron energy E at time instants ti and 
t2 marked on Fig. |6] for the bias corresponding to the NDR 
regime for (a) spin-down and (b) spin-up electrons. The posi- 
tion of the peak determines the energy of the spin- up (Eq^^, 

red line) and spin-down (Eq^^, blue line) resonance state lo- 
calized in the main QW. /je is the electrochemical potential 
of the emitter. 



electron interactions with the oscillating spin-down elec- 
tron density. 

The same effect, i.e., the coupling between the two 
spin-up resonance states, one localized in the emitter 
quantum well and the second localized in the QW, is 
responsible for the oscillations of the spin-up polarized 
current that occur in the NRD regime for this current 
component [cf. Fig. [2|a)]. This coupling leads to the os- 
cillatory exchange of the spin-up electrons between the 
emitter and the QW, which changes the tunneling con- 
ditions in an oscillatory manner. 



C. Effect of exchange interaction 

The results presented in Sec. |III| have been obtained 
without the exchange interaction [cf. Eq. (|3|]. We have 
also studied the effect of exchange interaction on the os- 
cillations of spin polarized currents. For this purpose 



we have estimated the exchange energy in the QW us- 
ing the formulsP^ = — (Sncr/Tr)^/^/^, which leads to 
■^ex ^ _Q meV for the time instant, for which charge 
density Ua- reaches the maximal value during the oscil- 
lation cycle. Inserting this formula to Eq. ([3| we have 
calculated the spin-dependent current- voltage character- 
istics and obtained very similar results to those presented 
in Fig. [2|a). In particular, we have found two bias 
regimes, in which THz current oscillations appear. The 
estimated frequency of these oscillations is smaller than 
that calculated without the exchange interaction, but 
still remains in the THz range. We conclude that the 
exchange interaction can lead to a slight modification of 
the current- volt age characteristics, but does not affect 
the predicted occurrence of the high-frequency spin cur- 
rent oscillations. 



V. CONCLUSIONS AND SUMMARY 

In the present paper, we have shown that - in the para- 
magnetic RTD under constant bias and external mag- 
netic field- the stable oscillations of the spin-polarized 
current can occur in the NRD regimes for the spin-up 
and spin-down current components. The analysis per- 
formed in Sec. [iVl reveals that these oscillations result 
from the coupling between the spin-up (spin-down) reso- 
nance state localized in the emitter region with the spin- 
up (spin-down) resonance state localized in the quantum 
well. This coupling leads to the oscillatory changes of 
both the electron density and the current in the nanos- 
tructure. 

In the RTD under the fixed bias and external mag- 
netic field, we deal with the following two types of the 
spin-polarized currents: (i) stationary (constant in time) 
current with the well-defined spin polarization, (ii) oscil- 
lating current with the oscillating spin polarization. We 
have demonstrated that at a certain bias in the NDR 
regime, at which the spin polarization of the net current 
reaches the maximum, the spin current polarization os- 
cillates with the THz frequency. 

Until now the THz oscillations of the current have 
been observed only in non-magnetic resonant tunnel- 
ing diodes^^^^ but they should be also occur in para- 
magnetic resonant tunneling diodes. We hope that the 
present results will stimulate the experimental search 
for this phenomena in the magnetic resonant tunneling 
structures. We would like to mention that the oscilla- 
tions of the spin-polarized currents are fully manifested 
only if the scattering on phonons and impurities is neg- 
ligibly small, i.e., this is the low-temperature effect. The 
scattering processes will damp the oscillations and finally 
lead to the steady current. 

The oscillations of the spin-polarized currents, pre- 
dicted in the present paper, can have important impli- 
cations for spintronic devices based on the paramagnetic 
RTD. The separation of both the spin components of the 
current is a basis for a fabrication of an effective spin filter 
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(selector), in which the spin polarization of the current 
is controlled by the bias voltage. However, if the para- 
magnetic RTD is designed to operate as a spin filter, we 
have to avoid the bias regimes, in which the current os- 
cillations appear, in order to obtain the well-defined spin 
polarization of the current. On the other hand, the cur- 
rent oscillations obtained under the constant bias can be 
applied to generate the spin-polarized current oscillating 
with the THz frequency. 

In summary, we have demonstrated that the spin- 
polarized current flowing through the paramagnetic RTD 
under the constant bias fixed in the NDR regime and 
in the constant external magnetic field can exhibit the 
THz oscillations. The oscillatory behavior of the spin- 
polarized current leads to important consequences for the 
operation of spintronic devices. On the one hand, this 



phenomenon is disadvantageous for the operation of spin 
filters as a parasitic effect that spoils the operation of spin 
filters. On the other hand, it can be used to design the 
generator of the spin-polarized currents oscillating with 
the THz frequency. 
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